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Implementations for quantum computing require fast single- and multi-qubit quantum gate op- 
erations. In the case of optically controlled quantum dot qubits theoretical designs for long-range 
two- or multi-qubit operations satisfying all the requirements in quantum computing are not yet 
available. We have developed a design for a fast, long-range two-qubit gate mediated by a photonic 
microcavity mode using excited states of the quantum dot-cavity system that addresses these needs. 
This design does not require identical qubits, it is compatible with available optically induced single 
qubit operations, and it advances opportunities for scalable architectures. We show that the gate 
fidelity can exceed 90% in experimentally accessible systems. 



I. INTRODUCTION 

Quantum information processing involves the manip- 
ulation of entanglement carried out by unitary gate op- 
erations between different quantum bits (qubits). Re- 
alistic quantum computing architectures require entan- 
gling gates between distant qubits. Optical photons pro- 
vide a natural vehicle to implement such interactions in 
many physical systems.^ As a result, architectures based 
on optically active qubits that can couple to photonic 
modes in optical cavities and waveguides, such as quan- 
tum dots, NV centers, and trapped ions are attractive 
for large scale quantum computing.'^ Quantum dots 
(QDs) in particular hold promise as qubits for such ar- 
chitectures, in part owing to their large dipole moments, 
which allow them to couple efficiently to the optical cav- 
ity modes and to photonic flying qubits for extended ar- 
chitectures. Qubits encoded by the spin of an electron 
in a QD have long coherence times which are five to six 
orders of magnitude longer than the typical picosecond 
scale of optical control. Successful initialization and read- 
out, as well as fast optical single s pin rotations, have 
been demonstrated in these systems.!^ In addition, im- 
portant advances have recently been achieved in work on 
coupled cavity-QD systems, including demonstrations of 
strong coupling and tunability, !^^'^^ and very recently 
full singie-qubit controL^i^ 

A critical step needed to advance the field is the de- 
sign of a two-qubit controlled gate operation mediated 
by an optical cavity mode. A viable two-qubit quantum 
gate requires that several criteria are met: (i) a long- 
range switchable physical interaction between qubits is 
available; (ii) the gate performs a unitary operation on 
one qubit depending on the state of the other qubit to 
provide a controlled operation; (iii) the operations are 
sufficiently fast compared to decoherence rates; (iv) the 
gate is compatible with single-qubit rotations (to form a 
universal set of gates); (v) the gate design is consistent 
with a multi-qubit system for scalability. 

So far, only local control of entanglement in closely 
spaced quantum dots (QD 'molecules') has been demon- 
strated experimentally.^^ For an experimental demon- 
stration of cavity-mediated entangling gates, a theoreti- 
cal design is needed that satisfies the above criteria, (i)- 



(v), while being experimentally simple and compatible 
with current technology. Existing proposals for cavity- 
mediated gates have not met these requirements; they 
are either incompatible with single-qubit gat es,'^^ limited 
to nearest-neighbor qubits, and/or require adiabaticity, 
either through adiabatic evolutiorP^or through adiabatic 
elimination of the auxiliary state.^^ As a result, they are 
much slower than what is needed from a quantum infor- 
mation processing perspective. Moreover, a careful as- 
sessment of the performance of such gates as a function 
of system parameters has not been given in the litera- 
ture, despite the key role it would play in experimental 
demonstrations . 

In this paper we give a novel design for an entan- 
gling control-z (CZ) two-qubit gatd^that satisfies all the 
above criteria. Our design does not require the QD ener- 
gies to be equal or dynamically tunable. As a result, our 
approach is compatible with single qubit operations and 
has a potential for many-qubit scalable architectures. We 
obtain fidelities in excess of 90% for realistic parameters. 
In the following we explain the concept of this all-optical 
gate, formulate the model, calculate the QD-cavity sys- 
tem spectrum, and analyze our design of the two-qubit 
gate protocol. The fidelity of the gate operation as a 
function of the system parameters is also calculated and 
provides a guide for experiment. 

II. TWO-QUBIT GATES 

The control-z gate is a maximally entangling two-qubit 
gate, and it is given by Ucz = diag(l, 1, 1, — 1). It is 
equivalent to the more familiar control-NOT (CNOT) op- 
eration up to single-qubit gates. Specifically, Ucnot — 
(l(8)H)J7cz(l«)H), where iJ = (J„i)/y2is the Hadamard 
gate. To see the entangling capability of CZ we can look 
at its action on a product state of two qubits. In partic- 
ular, when each qubit is in an equal superposition of the 
basis states, we have 

Ucz{\l) + |0)) ® (|1) + |0)) = 111) + |10> + |01> - |00), 

which is a maximally entangled two-qubit state, also 
known as a two-qubit 'cluster state'. Such a state is 
equivalent to a Bell state up to single-qubit rotations. 
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To implement the CZ gate, we need to accumulate a 
phase factor of —1 selectively to only one of the two-qubit 
basis states, taken to be |00) above. Meanwhile, to be 
able to perform single-qubit gates, the transition involv- 
ing state 1 00) and an auxiliary state should be performed 
in parallel with that involving |01) (or |10) for rotations of 
the second qubit) and its corresponding auxiliary state. 
To avoid dynamically tuning energies-a process that is 
costly in time and can compete with coherence times- 
we will use different classes of auxiliary states for single- 
qubit and two-qubit operations. In particular, we will 
use a near-resonance between the two-photon state of 
the cavity and the state where both QDs are excited. 

III. QUANTUM DOTS IN A CAVITY 

We focus on a system of two (singly) charged self- 
assembled InAs QDs in a photonic crystal microcav- 
ity. This structure can support in- and out-of-plane 
polarizations.!^ Due to strain the optical dipole transi- 
tion matrix elements in the InAs dots are anisotropic, 
resulting in efficient absorption of light with electric field 
polarization perpendicular to the QD growth axis. As 
a result, only the mode with electric field polarized in 
the plane of the crystal can be coupled to transitions in 
QDs. We take an external magnetic field to be applied 
in-plane (Voigt configuration), perpendicular to the QD 
growth direction. This will enable full single qubit con- 
trol as explained in Ref . fIU[ 

The system can be represented by two separate four- 
state QDs interacting with a single photon mode, as 
shown in Fig. [T|a). The two lowest energy states of 
each four-state system are the spin states of the elec- 
tron in each dot, which represent the qubit, |t) = c]^ ^|)o 
and ID — c|j^|)o, where n = 1,2 refers to the two 
dots and cj^^^ creates an electron of spin ^ (l) rela- 
tive to the uncharged QD state |)o. The two excited 
states in each dot are electron-exciton bound states, 
called trions (or charged excitons). They are complexes 
having total angular momentum 3/2. The two ±3/2 
states ('heavy holes') are energetically lower than the 
±1/2 ('light hole') states and thus form a pseudo spin 
Ifr) = 4,^|)o = ^.t^i.i'^LtDo ^^'^ similarly for 
where is the creation operator for a heavy hole. The 
trion states carry the pseudospin of the hole because the 
two electrons are in a spin singlet. We choose the spin 
quantization axis along the external magnetic field. 

The cavity couples to the trion transitions and pre- 
serves the (pseudo) spin orientation, |t) O |lt) and ||) O 
l-IJ.). In the rotating- wave approximation the cavity-dot 
interaction is 

Hqd-c = 5 X! {fn,-\'^n,\a + t\^^Cn,ia + h.c)j (1) 

n=l,2 

where a annihilates a photon in the cavity and g is the 
coupling between the trion transition and the cavity. We 




FIG. 1: The cavity-dot system, (a) Energies and relevant 
states of two QDs and cavity, (b-e) Interacting cavity-dot 
spectrum as a function of the cavity mode frequency, ujq. 

(b) Structure of crossings and corresponding states. Panels 

(c) and (d) show the anti-crossing splittings in the two- and 
single-excitation subspaces respectively. Panel (e) shows the 
energy structure of the qubit subspace, which is unaffected 
by the coupling to the cavity mode. The numbers in (c-e) 
give the states of the diagonalized Hamiltonian, and the f, i 
show the predominant spin character of each state far (to the 
right) of the avoided crossings. Vertical dashed lines indicate 
(ei +e2)/2. 



choose these coupling constants to be the same for the 
two QDs to simplify the presentation. This assumption 
is not important to the proposed procedure and can be 
relaxed when necessary. 

The spectrum of the cavity-QD system is shown in 
Fig. [T] as a function of the cavity frequency Wq- This 
representation does not suggest the need to tune luq dy- 
namically, but it helps to identify the region of optimal 
values. The spectrum is obtained by diagonalizing — 
Hqd + He + Hqd-c, where Hqd =J2n,i^e^i04i,i^ri,^ 

En.itld^n+uJhemtn,^, He = ujoa^ a, ^ =t,;, ^^(t) = 0, 
and 9{i) = 1. 

The Hamiltonian Hq conserves the total number of ex- 
citations. As a result the Hilbert space of the system 
separates into subspaces with different numbers of ex- 
citations. Each subspace contains several states, corre- 
sponding to different spin projections. The lowest set 
of four states (three energy levels) defines the two-qubit 
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subspace, |tt): lti)j lit); and HI), and has zero cav- 
ity photons; we call this the 'zero excitation' subspace. 
The other relevant subspaces are the 'one-excitation' sub- 
space, that has states with either one cavity photon or 
one trion, and the 'two-excitation' subspace, that has 
states with two trions (one per dot), states with one trion 
and one cavity photon, and states with two cavity pho- 
tons; see Appendix States in the 'one-excitation' part 
of the spectrum are approximately local to each quan- 
tum dot and interact with each other only very weakly, 
~ (g/Ae)^. They are the states that can be used for 
single-qubit control. The two-excitation subspace in- 
volves hybridized states of the two QDs and are ideal for 
a two-qubit gate. These states however are not directly 
accessible from the qubit subspace with a single pulse, so 
we make use of a series of control pulses. 

The laser pulses have momentum perpendicular to the 
photonic crystal plane to avoid Bragg shielding due to 
the photonic crystal. For definiteness we choose pulses 
with the same linear polarization as the cavity modej'^^l 

Vit) ^^np{t - tp)2 cos ujpt {Mnm\n){m\ + h.c.) . (2) 

The total Hamiltonian becomes 'H(t) ~ Ho + V(i), 
where Ho = HqU = I]„ (^^| and M„,„ = 

J2j=i,2,^{^\U^i^l,i^j,i + <^],i^^3,iW\^)- The subscript p 
enumerates the pulses used to perform the gate where 
each has frequency ujp and is centered at time tp. 



IV. IMPLEMENTATION OF CZ GATE 

The CZ gate has a simple diagonal form, which allows 
for a relatively straightforward design based on phases 
induced by resonant cyclic excitation of an auxiliary ex- 
cited state. The idea is to use the property of quantum 
two-level systems, in which a cyclic evolution from the 
ground state to the excited state and back to the ground 
state induces a minus sign to the latter. In the presence 
of additional, uncoupled states the minus sign is relative 
and thus constitutes a nontrivial quantum evolution. The 
pulse performing such an evolution is known as a '27r' 
pulse. Optical 27r pulses were proposed theoretically for 
single-qubit rotations in quantum dotd^ and two-qubit 
gates in quantum dot molecu lej^ and later used in their 
experimental demonstrations.'^E^ 

In our approach, the phase accumulation will be on 
state Iti), while keeping the phases of other basis states 
unchanged. This can be done by the following pulse se- 
quence: (i) a population inversion tt pulse, pulse A, tuned 
to transition cji = loa = E^ — Eq between qubit state |tt) 
and the excited state with similar spin configuration, see 
Fig[IJd)-(e). The pulse is also in resonance with Eq — E2 
transition, and thus it creates a trion in the first QD only: 
both Itt) and |ti) are transformed in the same way and 
accumulate a phase factor of —i each, (ii) A 27r, or phase, 
pulse (pulse B) with frequency UI2 = lob = Eiq — E4, 



see Fig [Tjc)-(d). This induces a transition between the 
'one-excitation' states previously created and one of the 
'two-excitation' states. Note that if g = or we are far 
detuned, w ^ Ae, the transition Eiq — E2 would also 
occur. This would correspond to a single qubit opera- 
tion on the second qubit, i.e., Ifft) and would both 
acquire a phase factor of —1. A nonzero g induces forma- 
tion of two-excitation states that have different energies, 
c.f. the energy of state 16 and the sum of energies of 
states 4 and 8. As a result, the state fft: or state 4, ac- 
quires the factor of —1 after the pulse, while state ||t) 
does not. (iii) Finally, we apply the population inversion 
pulse A again, ^3 = uja, to restore the system to the 
qubit subspace. This gives additional factors of —i to 
both Itt) and |tt), as mentioned above. The two phase 
factors of (— z) induce a minus sign to states |tt) and |tt)j 
while the 2tt pulse cancels that sign in state |tt)- The 
phase between the control pulses A and B does not enter 
the result and therefore pulses with unequal frequencies 
do not have to be phase locked, which is a significant 
experimental convenience. 

A physical explanation of this approach is the follow- 
ing: because each QD is off-resonant from the cavity, 
when only one of the QDs is excited and no other exci- 
tations are present in the system the excited QD can be 
roughly thought of as isolated, i.e., decoupled from the 
cavity and from the other QD. Thus, single excitations 
can implement single-qubit operations without disturb- 
ing the rest of the system. On the other hand, when both 
QDs are excited they are closer to the resonance with the 
cavity state. As a result, there is a large mixing between 
cavity states and the states of both QDs. Thus, using the 
two-excitation regime is a natural venue for performing 
two-qubit conditional operations while maintaining the 
ability to manipulate each QD spin separately. 



V. FIDELITY 

Now we consider the gate fidelity, which is a measure 
of how close our operation is to the target gate. There 
are two types of fidelity losses, those caused by unin- 
tended coherent dynamics due to coupling of the lasers 
to off-resonance transitions and those originating from 
random processes such as trion recombination. First we 
focus on the former mechanism. The unintended tran- 
sitions can cause Ug to deviate from the ideal Ucz and 
effectively cause loss of coherence in the qubit-subspace, 
even though the entire operation involving excited states 
is unitary and coherent. 

To analyze this type of decoherence we compute the 
average fidelity, F, of the gate operation, as explained in 
detail in Appendix [B| including transitions 0-4, 4-6, 2-5, 
3-7, 0-8, 2-10, 1-9, 3-11 for pulse A and transitions 4-16, 
6-18, 0-8, 2-10, 1-9, 3-11 for pulse B. Other transitions 
are negligible either due to vanishing matrix elements 
or to large detuning. We chose different pulse widths 
for pulse A and pulse B aA = 2a b = 2a. In Fig. [2](a) 
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FIG. 2: Fidelity of CZ gate with imperfections resulting from coupling of the pulses to neighboring off-resonance optical 
transitions (a) as a function of ljq for a/oje = 0.1 for Ae/uJe = 8.33, 16.67, 25.00, 33.33, as indicated, and (b) as a function of 
the spectral separation Ae between the QDs for different values of the pulse bandwidth, a/u^ = 0.01, 0.02, 0.1, 0.15, 0.3 
as indicated by the dashed arrow. Each point is computed for the optimal value of loq from Fig. [2][a). The vertical lines mark 
the values of 0J2 from panel (a). In both panels (a) and (b) we used g/uJe = 3.33 and ujh = tUe/3. 




FIG. 3: (a) Fidelity of the two-qubit CZ gate in presence of decoherence due to trion recombination and cavity decay. The 
fidelity is plotted as the function of the trion decay time and the cavity mode quality factor, (b) The temporal profiles of the 
pulse sequence for tj/oje = 0.2, luh ~ <^e/3, Ae/oje = 8.33, g/oJe ~ 3.33, and cje ~ 0.12meV. 



the fidelity is plotted as a function of the difference be- 
tween the cavity mode frequency and the transition 
frequency of QDl ei for varying values of the frequency 
of QD2 £2. The qualitative features of the plots can be 
understood as follows: when the cavity mode frequency 
is much smaller or much larger than the QD frequencies, 
QD-cavity hybridization is negligible, and we are in a 
regime of two independent qubits. This causes attenua- 
tion of fidelity towards both sides of the plot. The dip 
in the middle occurs because, as the cavity is tuned, the 
target transition of pulse B (transition 4-16) becomes de- 
generate with transition 3-11, and therefore state | 4,4,) is 
also affected by pulse B, resulting in strong unintended 
dynamics. Note that at its high values the fidelity does 
not vary strongly with ei and £2- As a result, gates be- 
tween several different pairs of quantum dot spin qubits 
can be performed with high fidelity using only one cavity 
mode to mediate the interactions, which is an intriguing 
opportunity for scalable architectures. 

Fig. ^h) shows the fidelity as a function of the spectral 
separation A£ between the trions in QDl and in QD2 for 
different pulse bandwidths a. When A£ is small (com- 
parable to We) the fidelity drops appreciably. This drop 
is the result of coupling in the 'one-excitation' subspace, 



i.e., the assumption that an excited QD is isolated from 
the rest of the system is no longer valid. Thus it also 
identifies the regime where fast optical single-qubit con- 
trol is not possible. In the region of higher fidelities, 
where Ae/cje > 10, the fidelity approaches its maximal 
value for longer pulses and starts decreasing more rapidly 
for a/uje ^ 0.2 due to involvement of a larger number of 
unintended transitions. 

Next, we consider the effects of decoherence due to 
losses during the gate. The main contributions come 
from trion recombination and cavity photon leakage. The 
typical linewidth of the trion state, Ftr, in InAs QDs is 
~ IfieVp^ The loss rate associated with the cavity is 
Tc = loq/Q- State-of-the-art microcavitie^l^can have Q's 
up to ^ 10^, which gives Tc ~ 10/xeV. We calculate t he fi- 
delity using the standard master equation formalisrrPIIlIl 
and include states from to 19, see Appendix [BJ The 
fidelity as a function of Q and l/Ftr in shown in Fig. 3] 
It is maximized when the pulses overlap to reduce the 
time the excited states are occupied. We choose luq from 
the maximal fidelities, as in Fig. [2](a) for each point of 
Fig.|3] We see that fidelities in excess of 90% are possible 
for realistic values of the parameters. 
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VI. CONCLUSIONS 

In summary, we have developed a design for a cavity- 
mediated entangling gate between two spin qubits that 
satisfies the criteria for a realistic two-qubit operation. 
Our control-z gate is compatible with available single 
qubit operations and with natural inhomogeneities in op- 
tical resonances. It can thus accommodate several qubits 
that couple pairwise with appropriate control laser fre- 
quencies, opening a path to scalable architectures. It 
may also be useful for hybrid quantum computing ap- 
proaches with various physical systems.^^SlWe have shown 
that the gate fidelity is at least 90% for current experi- 
mental parameters. Higher fidelities can be achiev ed in 
various ways such as using pulse shaping technique^^HMI 
and engineering higher finesse cavities. 
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Appendix A: The spectrum 

In the rotating wave approximation the Hamiltonian 
[Eq. (1) from the main text] conserves the total number 
of excitations. Therefore it can be diagonalized indepen- 
dently in each excitation-number subspace. The lowest 
energy set of four states (three energy levels) corresponds 
to a subset with zero excitations. It represents the two- 
qubit subspace with zero cavity photons, 

0^ in>|o),i^ in)|o),2^ |;t)|0),3^ |U)|0), (Al) 

where |0) is the vacuum state of the cavity. The corre- 
sponding energies are controlled by the magnetic field via 
Zeeman splitting. For typical values of magnetic field of 
~ 1 T used in the initialization and readout and single- 
qubit experiments the splitting between Eq,E^ and E12 
is ^ 0.1 meV. The micro-cavity optical mode is coupled 
to the excitonic transitions in each quantum dot with the 
transition energies ~ eV. As a result, the qubit subspace 
is not affected by the cavity. 

The one-excitation subspace occurs at the optical fre- 
quency, ~ eV from the qubit subspace energies: 



dot 1 
dot 2 
cavity 



|^t)|0),|^t)|0),|M)|0),|4i)|0), (A2) 
|n)|0),|n)|0),|W|0),|U)|0), (A3) 

itt)|i),in>|i>,i;t>|i>,iu>|i>, (A4) 



where |1) denotes the state with a single photon in the 
cavity. The energy gap Ae between states ( A2 ) and ( A3 ) 



is due to the fact that the two dots are not identical in 
size, shape, and strain environment, which affects the 
excitonic transitions. The typical variation in trion tran- 
sition energies is ~ 1 — 20 meV. The energy of the cavity 



mode. Wo, is fixed during the gate operation but can be 
set to an optimal value during sample growth. The in- 
teraction with a cavity photon shifts the energies and 
mixes trion and photon states. The energies of the re- 
sulting states can be found analytically: note that states 
( A2|| A4 ) are always coupled in triplets. For example, the 
state7tt)|l) interacts only with |trt)|0) and |tfr)|0). For 
each triplet we have 

where ^ =t or i, e„_^ = e„ and = e„ + iVh - c^e- 
Each triplet forms two anti-crossings when ljq is swapped 
across the trion energies [see Fig. 1(b) and Fig. 1(d) of 
the main text]. When g ~ Ae or g ^ Ae, the two excited 
quantum dot states can mix and form spin-entangled 
states. For experimentally accessible systems of quantum 
dots in a micro-cavity the coupling strength g is substan- 
tially smaller than the variation in trion energies Ae and 
the mixing is negligible. In the limit g <^ Ae the interac- 
tion between one-excitation states from different QDs can 
be estimated by analyzing the difference 6uj^ in transition 
energies between : |tt) IfTt) and uj'^ : |ti) Ifti). 
From Eq. (A5| we find uj^ = e + g"^ / /{uja, Ae) and 



X — UJQ 



e + gyfi^A^^^ + We - ujh), where fiy,x) 
g'^/{x — e + y). Since 



cj/i ~ g ^ Ae it 

is easy to show that Soj^ ^ —g'^uJe/ Ae'^ . This should be 
compared to the typical inverse lifetime of the trion state, 
~ IfieV (in energy units) or ~ We/lOO. As a result for 
LUe/ Ae ~ 10, uj^ and are practically indistinguishable. 
This result is confirmed numerically by computing the 
spectrum (and the states) for different values of Ae. It 
also holds for other transitions between the qubit and the 
one-excitation subspace states. Therefore we conclude 
that the one excitation subspace cannot be used for a 
two-qubit operations. It can, however, be used to per- 
form fast single qubit operations as described in Ref. [23] 
by using the localized trion state. 

In order to find useful non-local states that can mediate 
a two-qubit gate we investigate the two-excitation sub- 
space. In this subspace the states are coupled in groups 
of four, e.g. |^^,0), |^t,l), in,l), in,2): 

{e2,i+ei^^-E){uj+ei^^~E){ujo+e2.i-E){2u}-E) (A6) 
^g\e2,^+ei,^ + 2uQ^2Ef. 

The spectrum has a more complex structure, see Fig. 1(c) 
of the main text. The two-excitations subspace provides 
non-local quantum-dot-cavity states, such as state 16, 
which has two trions (one in each dot). The energy of 
such state is different from the combined energy of two 
trion states localized in each dot, such as E4 and Eg,, 



AEie^4 7^ AE^fl + AEs 



(A7) 



where AEn,n' = En ~ En' ■ This is the basis for the 
two-qubit conditional phase gate in this work. Us- 
ing a perturbative approach like that above, we obtain 
A£;i6.4 - {AEifl + AEs,o) ~ -g^Ae. 
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Appendix B: Gate Fidelity 

The fidelity of the gate described in the main text is 
affected by two type of processes: (i) induced unintended 
transitions between the states of the qubit-cavity system 
and (ii) real losses due to cavity leakage and trion re- 
combination. We first estimate losses due to unintended 
but coherent dynamics. We include transitions 0-4, 4-6, 
2-5, 3-7, 0-8, 2-10, 1-9, 3-11 for pulse A, and 4-16, 6-18, 
0-8, 2-10, 1-9, 3-11 for pulse B (w^). Other transitions 
are negligible either due to vanishing matrix elements or 
to large detuning. We compute the wave function after 
the A-B-A pulse sequence for each basis configuration of 
the qubit subspace as initial state (evolution is linear and 
therefore the resultant wave function for any initial qubit 
state can be easily recovered). To simplify calculations 
here we resort to analytically solvable Rosen-Zener pulse 
shapes,!^ i.e. ^p{t) = rJpSech(f7pi) with a a = 2crs — 2f7, 
to calculated transition amplitudes and phases for reso- 
nant and off-resonance transitions for each pulse. Given 
the initial, and final, = U\il}o), wave function, 

the fidelity can be computed as 



(Bl) 



where U]j^ is the evolution operator corresponding to 
the ideal CZ gate. The value of F{ipQ, ip) depends on the 
initial state of the two-qubit system and therefore can 
vary depending on the choice of algorithm and initial 
data. We therefore compute the average fidelity F by 
taking average over all possible initial states of the two- 
qubit system. 



(B2) 



— — 



ijnm 



20 



n\UlzU\i){j\U^Ucz\rn} 



The integration J dtpQ is performed over all complex am- 
plitudes that define the initial state in the basis |«), and 
n, m run over all basis states ftj tij it; The re- 
sults are presented in Fig. 2 and the discussion is given 
in the main text. 

In order to account for both unintended dynamics and 
actual losses we have to calculate the reduced density 
matrix, p{t), of the two-qubit sub-system for the duration 
of the pulse sequence. The reduced density matrix can be 
found within the Bloch-Redfield master-equation (ME) 
formalism 



ip^ [H + Vit),p] 



(B3) 



p + p 



PlP. 



where Ps = |/s)(«s|, \is) and |/s) are initial and finial 
states (in the spin basis) corresponding to the s-th de- 
cay process with rate F^. Solving the above equation 
directly is computationally involving due to the presence 
of two different time scales: fast, associated with the 
laser driving frequency, and slow, coming from the time- 
dependence of the pulse shaped envelope. To simplify 
the computation we transform the ME to the eigenbasis 
of H and use the rotating wave approximation. 



(B4) 



VspVl - 



VlVsP + pVlVs 



where p = e^'^^W^ pUe-''^o\ Ps = PsU and 



Vit) 



V(t)e 



-iUot 



Note that the trion decay pro- 



cesses can involve photons with any in-plane polariza- 
tion (along or perpendicular to the applied magnetic 
field). Therefore, for each trion we have Fs = Ftr: 
Ps ^ {|t)(^l,i;)(^|,|t)(4|,|;)(4|}. Leakage of pho- 
tons from the cavity is modeled as Tg = F^: Ps — >■ 
{|0) (1|, \/2|l)(2|, etc.}. Due to additional (pseudo)spin- 
fiip electron-hole recombination processes, more states 
are involved than for the coherent case discussed above 
and we include states from to 19. We chose to use Gaus- 
sian pulse shapes ilp(t) = {Vlp/ ^Jtt /2) exp{—2t'^ap/n'^} 
for numerical convenience and apply the same pulse se- 
quence as before with aA — ^.as — 2(7. 

Since a separable quantum wave function is no longer 
accessible, fidelity has to be defined differently. 



F{ijo,pM) = \/{i^o\ulzPUcz\i^o) 



(B5) 



In this case the average fidelity is computed as 

' + ^'^^''^{n\ul,p{\^){J\}UczH (B6) 



ijnm— {1,4:} 
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which is the generalization of Eq. (B2| for the case of 



non-unitary evolution of pure initial state. This is possi- 
ble due to the fact that the evolution of the density ma- 
trix is still described by a linear (but non-unitary) super- 
operator, i.e. p{t) = T exp{—i J^dtLnit) — tL)|V'o)(^o|j 
where L^O — [H,0]. As a result, the complex coef- 
ficients that define initial {\tpo)) and target (Uczli^o)) 
states in the basis \i) can be integrated out in exactly 
the same way as for Eq. ( B2 ) . The results are presented 



in Figs. 2 and 3 in the main text. 
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